Density Matrix Renormalisation Group Calculations for Two-Dimensional Lattices: 
An Application to the Spin-Half and Spin-One Square-Lattice Heisenberg Models 
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A new density matrix renormalisation group (DMRG) approach is presented for quantum systems 
of two spatial dimensions. In particular, it is shown that it is possible to create a multi-chain-type 
2D DMRG approach which utilises previously determined system and environment blocks at all 
points. One firstly builds up effective quasi-lD system and environment blocks of width L and these 
quasi-lD blocks are then used to as the initial building-blocks of a new 2D infinite-lattice algorithm. 
This algorithm is found to be competitive with those results of previous 2D DMRG algorithms and 
also of the best of other approximate methods. An illustration of this is given for the spin-half and 
spin-one Heisenberg models on the square lattice. The best results for the ground-state energies per 
bond of the spin-half and spin-one square-lattice Heisenberg antiferromagnets for the N — 20 x 20 
lattice using this treatment are given by E g /NB ~ —0.3321 and Eg/Ns = —1.1525, respectively. 

PACS numbers: 75.10.Jm, 75.50Ee, 03.65.Ca 



I. INTRODUCTION 

The subject of one-dimensional lattice quantum spin 
systems at zero-temperature contains a large number of 
exact solutions, for example, via the Bethe Ansatz so- 
lution 0, 0, 0, 0. However, such exact solutions have 
not as yet been universally obtained for systems of quan- 
tum spin number, s > 1/2, or for lattices of larger than 
one spatial dimension, or indeed for frustrated spin sys- 
tems. It is noted however that exact solutions for such 
systems do exist for various specialised cases. Recent 
density matrix renormalisation group (DMRG) calcula- 
tions |5j have been extremely useful in determining the 
ground- and excited-state properties of a whole host of 
one-dimensional (ID) or quasi-lD spin systems to ex- 
tremely high accuracy. The DMRG method is not lim- 
ited by the magnitude of the quantum spin number or 
by the presence of frustration, and DMRG has been used 
to determine the non-zero temperature and dynamical 
properties of many quantum systems. 

It is perhaps still fair to say that the full power of 
the DMRG method has largely been restricted to one- 
dimensional systems, although recent calculations have 
very successfully extended |1 Q, H H [H EJ 13 pre- 
vious DMRG calculations for one-dimensional quantum 
systems to two spatial dimensions (2D). Perhaps the sim- 
plest such 2D DMRG approach is the "multi-chain" ap- 
proach [f| 0,0 EI m which, during the infinite-lattice al- 
gorithm, ones keeps the width of the lattice constant and 
then one grows the lattice height by adding whole rows 
or partial rows of sites. The multi-chain finite-lattice al- 
gorithm may then be implemented and one repeatedly 
sweeps through all lattice sites until a fixed point of the 
renormalisation group (RG) process is reached. This is 
in direct analogy with the traditional one-dimensional al- 
gorithm, and so (in some sense) one has mapped the 2D 
problem onto an effective ID problem. This multi-chain 
approach has been seen to work well for relatively small- 



sized lattices in Ref. ^lj in which whole lines of sites 
are added at a time. Indeed, a "finite-lattice sweep" was 
performed after each addition of a line of sites (and not 
just after the final N = LxL lattice had been created) in 
this calculation , and this was seen to greatly enhance 
the accuracy of the method. 

Another approach 10j was to split the two-dimensional 
lattice into four pieces, namely, two-dimensional system 
and environment blocks separated by a lines of spins 
which itself was split into distinct blocks, namely, one- 
dimensional system and environment blocks, and free 
spin(s). This "four block" approach is a truly two- 
dimensional algorithm in the sense that one is able to 
grow the two-dimensional lattice in a controlled manner 
using the results (i.e., ID and 2D blocks) of previous iter- 
ations and a lattice of arbitrary size may be thus finally 
obtained. However, this approach suffers from a num- 
ber of problems such as, for example, the need to diago- 
nalise very large density matrices for the two-dimensional 
blocks for relatively modest numbers of states. We also 
note that these density matrices for the 2D blocks may 
also be highly singular. 

Another recent and very successful 2D DMRG ap- 
proach used two blocks at opposite corners of a square 
lattice with two "free" spins placed at particular posi- 
tions along the boundary between the system and envi- 
ronment blocks. It was seen for this method that there 
was no need to set the width of the lattice to be constant 
and that previously defined blocks for a given N = Lx L 
lattice could be used in the following DMRG iterations to 
form lattices of size N' = L+lxL + 1. Furthermore, this 
approach was found to be very efficient for the (relatively 
small) lattices (N = L x L with L < 12) considered, and 
excellent results for the ground-state energy per bond of 
the spin-half Heisenberg model on the square and tri- 
angular lattices were obtained. An application of the 
DMRG method for highly anisotropic two-dimensional 
systems was also recently performed |12|. 
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In this article a new 2D algorithm is presented which 
builds up two-dimensional lattices in a fast, simple and 
efficient manner. Quasi-f D blocks of up to a given length, 
L, are determined and thereafter the width is kept con- 
stant - in analogy with previous multi-chain approaches. 
However, these quasi-lD blocks are then used to form the 
first initial steps of an algorithm which uses 2D blocks of 
various shape and which have been determined during 
earlier DMRG steps in order to finally build up a lattice 
of size N = L x L. Lattices of arbitrary size may thus 
be considered, although this 'size' is pre-set at the start 
of the DMRG calculation. This is in marked contrast to 
previous multi-chain approaches. The 2D "finite-lattice" 
DMRG algorithm has also been implemented here in or- 
der to arrive at a fixed point of the renormalisation group 
process, and this part of the algorithm is similar (if not 
identical) to previous multi-chain algorithms. 

Note that the details of the DMRG approach are de- 
ferred until the appendix. However, it is important to 
note that the current approach does not utilise a row of 
sites or partial row in order to build up the lattice in the 
infinite-lattice algorithm. Rather, it does use previously 
defined blocks for the environment block which have both 
width and height at all points. The total system size is 
thus incremented by more than a single line of sites at a 
time, although the width is kept constant. 
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FIG. 1: The DMRG ground-state energies per bond of the 
spin-half Heisenberg on a JV = 6x6 square lattice for 
M = {50, 100, 150}. Note that each "iteration" indicates a 2D 
DMRG "move" described in detail in the appendix. Note that 
the infinite-lattice algorithm ends after 22 iterations (close to 
the peak in the Eg/Ns) and the results thereafter refer to 
those for the finite-lattice algorithm. 



II. THE HEISENBERG MODEL FOR THE 
SQUARE LATTICE 

The Heisenberg model is given by 

H = J2 S * ■ s o > (!) 

and where (i, j) runs over each nearest- neighbour bond 
on the square lattice (counting each bond once and once 
only) and open boundary conditions are assumed. The 
results of the best of other approximate methods for the 
spin-half case [H [H E []J predict that E g /N B « 
—0.3347 (where N B is the number of bonds) and that 
61-62% of the classical ordering remains for the quantum 
system in the infinite-lattice limit, N — > oo. The results 
of the best of other approximate methods for the spin-one 
case [H 13 predict that E g /N B » -1.164 (where N B 
is the number of bonds) and that 80-81% of the classical 
ordering remains for the quantum system in the infinite- 
lattice limit. 



III. RESULTS 

A plot of the ground-state energy per bond for the s = 
1/2 Heisenberg antiferromagnet for the N = 6 x 6 lattice 
for M = {50, 100, 150} is presented in Fig. 2] It is seen 
that the trend of the ground-state energy is to decrease 
with the number of iterations and the minimal solution is 



quickly reached. Note however that edge effects are seen 
where, for example, the ground-state energy is slightly 
lower for moves 1-3 of the algorithm presented above. 
Also, the infinite-lattice algorithm ends after 22 "moves" 
indicated on the abscissa as "iterations" in Fig. Q]and we 
note that the ground-state energy per bond has a peak 
at this point. Thus, the finite-lattice algorithm reduces 
the ground-state energy per bond by about 6% to 8%. 
Slightly larger gains (up to about 10%) from the finite- 
lattice algorithm are possible for the larger lattices. 

A value for the ground-state energy is taken from one 
of the sites near to the centre of the 2D square lattice 
for all of the N = L x L lattices considered here. By 
contrast, it is noted that the number of iterations needed 
to obtain convergence increases strongly with increasing 
N. It is observed that the number of iterations neces- 
sary can be as large as 20000 or more for L = 12 and 
L = 20. Indeed, it was seen that the ground-state en- 
ergy sometimes decreases in a step-like fashion (even af- 
ter 10000 or more iterations large step-like changes were 
sometimes observed). An example of this is given in a 
plot of the ground-state energy per bond for the s — 1 
Heisenberg antiferromagnet for the N — 12 x 12 lattice 
for M = {50, 100, 150} of Fig. It would be inter- 
esting to test whether such step-wise changes occur for 
other 2D DMRG algorithms and if this were a character- 
istic of implementing the DMRG method in two spatial 
dimensions. Finally, an explicitly assumption of these 
simulations was that they had been run long enough in 
order to reach the RG fixed-point. 

Note that the infinite-lattice algorithm ends after 



3 



-1.06 


1 i | 






-1.08 




— M=50 

s~\ s~\ \ it i AA 

O— O M=1UU 
x-x M=150 




-1.1 

C /IN- 
g B. L12 






- 


-1.14 








-1.16 






I.I.I. 



I I I I I I I 

5000 10000 15000 

Number of Iterations 



FIG. 2: The DMRG ground-state energies per bond of the 
spin-one Heisenberg on a N = 12 x 12 square lattice for 
M — {50, 100, 150}. Note that each "iteration" indicates a 
2D DMRG "move" described in detail in the appendix. Note 
that the infinite-lattice algorithm ends at the peak in E 3 /Nb 
and the results thereafter refer to those for the finite-lattice 
algorithm. 



about 80 "moves" in Fig. [3 (which are again indicated on 
the abscissa as "iterations" ) and we note that the ground- 
state energy per bond has a peak at this point. Thus, the 
finite-lattice algorithm reduces the ground-state energy 
per bond by about 8% to 10%. 

These results may also be compared to results of pre- 
vious DMRG calculations [13 in Tab. 0for the spin-half 
Heisenberg model on square lattices of varying size. Two 
algorithms were used in Ref. , which were, namely, a 
2D algorithm in which the 2D system and environment 
blocks are grown from the corners and a multi-chain ap- 
proach (also described in Refs. 0013) m which whole 
lines of sites were added at a time during the infinite- 
lattice approach. The results of these approaches shall 
henceforth be referred in this article to as DMRG (a) 
and DMRG (b) results, respectively, and the results of 
the present approach using the infinite-lattice algorithm 
shall be referred to as DMRG (c) results. 

It is seen from Tab. U]that we obtain reasonable cor- 
respondence between the present DMRG (c) results and 
the previous DMRG (a) and (b) results for all the lattices 
considered. A simple 'heuristic' extrapolation of these re- 
sults in the limit M — > oo by assuming that the ground- 
state energy scales as a "power-law" with respect to M, 
which ought to be the case for a gapless system. The 
best fit of such a power-law dependence was determined 
and the extrapolated values were determined. Note that 
this simple treatment gave an extrapolated value for the 
ground-state energy per bond of E g /N B = —0.36219 for 
the N = 6 x 6 lattice. This result is in good agreement 
with the extrapolated value for the ground-state energy 
of the spin-half square-lattice Heisenberg antiferromag- 
net of Ref. [TJ, namely, E g /N B = -0.36212. The ex- 



TABLE I: DMRG results for the ground-state energy per 
bond of the spin-half Heisenberg antiferromagnet on the 
square lattice for lattices of size N = L x L with 
L = {6, 12, 20} and for number of DMRG states M = 
{50, 100, 150}. The new results are referred to as DMRG 
(c) results and these results are compared to two previous 
DMRG two-dimensional approaches (see Ref. [Tl|') referred 
to as DMRG (a) and DMRG (b) results, respectively. It is 
seen that the present DMRG results are in reasonable agree- 
ment with these previous 2D DMRG results - especially for 
the N = 12 x 12 lattice. 





N 


M = 50 


M = 100 


M = 150 


DMRG (a) 


6x6 


-0.361972 


-0.362096 




DMRG (b) 


6x6 


-0.361919 


-0.362089 




DMRG (c) 


6x6 


-0.35847 


-0.36135 


-0.36184 


DMRG (a) 


12 x 12 


-0.337374 


-0.341588 




DMRG (b) 


12 x 12 


-0.332574 


-0.338833 




DMRG (c) 


12 x 12 


-0.33260 


-0.33757 


-0.33984 


DMRG (c) 


20 x 20 


-0.31763 


-0.32441 


-0.32679 



trapolated result for the ground-state energy per bond of 
the spin-half antiferromagnet on the N = 20 x 20 lattice 
is given by E g /N B = -0.3321. This result for the N = 
20 x 20 lattice is in reasonable agreement with the results 
of the best of other approximate methods [l3L HEL Efl . 
These methods predict that E g /N B ~ -0.3347 in the 
infinite-lattice limit, N — > oo, and note that this result 
ought to be quite close to the result of the N = 20 x 20 
lattice, even with closed boundary conditions. (Note that 
no extrapolations of the 2D DMRG results in the limit 
./V — > oo were determined due to the small number of 
data points and the fact that a relatively large lattice 
was considered anyway, although such an extrapolation 
is possible to perform.) 

Results for the spin-one Heisenberg model are also 
shown in Tab. [H]for square lattices of size N = 12 x 12 
and N = 20 x 20. These these results were again ex- 
trapolated in the limit M — > oo and we obtained an 
extrapolated result of E g /N B = —1.1525 for the spin- 
one antiferromagnet on the N = 20 x 20 lattice. This 
result for TV = 20 x 20 the lattice is seen to be in mod- 
est agreement with those results of the best of other ap- 
proximate methods 0, for both the ground-state en- 
ergy per bond, evaluated in the limit N — > oo, namely, 
E g /N B » -1.164. 

It may also be seen from Tab. [I] that the DMRG 
(a) results of the 2D algorithm of Ref. 11] appear to 
perform better in terms of the ground-state energy per 
bond of the spin-half Heisenberg model than either the 
DMRG (b) and the present DMRG (c) approaches for 
all of the lattices considered. It is thus interesting to 
consider why different DMRG approaches for exactly the 
same model and lattice at equivalent levels of approxima- 
tion yield differing results. It is known from the sim- 
ulation of one dimensional chains with periodic boundary 
conditions that the ground state does not become fully 
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TABLE II: New DMRG results for the ground-state energy 
per bond of the spin-one Heisenberg antiferromagnet on the 
square lattice for lattices of size N = L x L with L = {12, 20} 
and for number of DMRG states M = {50, 100, 150}. The 
DMRG results are compared to those results of SWT JT^ . 
cumulant series expansions [l5|. and the CCM for the 
infinite-lattice limit. 





N 


E a /N B 


DMRG(c) M = 50 


12 x 12 


-1.15472 


DMRG(c) M = 100 


12 x 12 


-1.16780 


DMRG(c) M = 150 


12 x 12 


-1.1715 


DMRG(c) M = 50 


20 x 20 


-1.10654 


DMRG(c) M = 100 


20 x 20 


-1.13939 


DMRG(c) M = 150 


20 x 20 


-1.1462 


CCM 


oo 


-1.1646 


SWT 


oo 


-1.1641 


Series Expansions 


oo 


-1.16395(1) 



were found to be a reasonable agreement with both previ- 
ous DMRG results for the spin-half square-lattice model 
[l"ll | with M = 50 and M — 100. Furthermore, the new 
DMRG results for the N — 20 x 20 lattice were seen to 
be in reasonable agreement with the results of the best 
of other approximate methods 0, 0, 0, 0, 0] de- 
termined in the limit N — > oo. These new results, and 
those results of previous and highly successful DMRG 
calculations in 2D @, 0, H, El E3 > now snow that the 
DMRG method is quickly becoming a powerful tool in 
order to simulate and understand strongly interacting 
two-dimensional quantum many-body systems. 



translationally invariant, thus indicating that the "global 
minimum" obtained may be somewhat dependent on the 
build up considered. However, even for the same topo- 
logical buildup and the same final number of states kept, 
the outcome may depend qualitatively on which way the 
number of states kept was increased during the sweeping 
process of the finite-lattice algorithm J2jj • 



IV. CONCLUSIONS 

In this article a new density matrix renormalisa- 
tion group (DMRG) algorithm was presented which 
determines the properties of two-dimensional square- 
lattice spin-half and spin-one Heisenberg antiferromag- 
nets. This algorithm builds up 2D lattices of arbitrary 
size in a straightforward manner, and, although concep- 
tually similar to earlier multi-chain approaches for 2D 
lattices, the method differs from these earlier approaches 
in practice because one firstly builds up effective quasi- 
1D system and environment blocks of width L. These 
quasi- ID blocks were then used to form the initial steps 
of a infinite-lattice algorithm in order to build a lattice 
of arbitrary, although (as for the ID DMRG algorithm) 
preset at the start of the calculation, size. It has thus 
been proven that it is possible to construct a multi-chain 
DMRG approach for 2D lattices which uses previously 
determined system and environment blocks at all points. 

This new 2D DMRG approach tested for the spin-half 
and spin-one square-lattice Heisenberg models in order 
to predict their ground-state energies with great success. 
The best results for ground-state energy of the square- 
lattice Heisenberg model for the N = 20 x 20 lattice for 
the spin- half and spin-one models were found E q /Nb = 
-0.3321 and E g /N B = -1.1525. The results for the 
Heisenberg model on the N = 6 x 6 and N = 12 x 12 



Future 2D DMRG calculations for the spin-half and 
spin-one Heisenberg models are envisaged for even larger 
square lattices and for larger values for M. Other future 
calculations using the new DMRG approach could also be 
for the spin-half and spin-one Heisenberg models on the 
triangular and Kagome lattices. Indeed, extensions of the 
approach presented in this article to lattice quantum spin 
systems with frustrating next-nearest-neighbour bonds 
are also possible. It is finally noted that the 2D algo- 
rithms considered in this article may be applied to a wide 
range of both fermionic and bosonic systems. 
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in performing DMRG calculations. 
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FIG. 3: The one-dimensional system block is the lower block 
of sites and it is used to initialise the infinite-lattice algorithm. 
The "free" spins SI and S2 are indicated by the unfilled cir- 
cles and the environment is the upper block (chain) of sites. 
Diagram (a) thus shows a one-dimensional system block (6, 1) 
of width 1\d(~ 6) and of height 1. Figure (b) shows the sys- 
tem block (7, 2) which is of width L(— 9) and height 2. (Note 
that the differently shaped blocks are uniquely described by 
the indices {x, y) shown above.) 



APPENDIX A: DENSITY MATRIX 
RENORMALISATION GROUP (DMRG) 
ALGORITHM FOR TWO SPATIAL DIMENSIONS 



The method of building up the 2D lattice is now in- 
troduced and the first step in this algorithm is to build 
up one- dimensional lattices of width L and height of one 
and two sites (shown in Fig. |3J). These blocks are used as 
the initial starting point for the two-dimensional infinite- 
lattice algorithm. Note that each differently shaped block 
in Fig. 131 is uniquely described by the indices {x,y), as 
shown in the figure and that this labelling of the dif- 
ferently shaped blocks is crucial to this 2D implementa- 
tion of the DMRG algorithm. System blocks are thus 
uniquely denoted by the indices x and y, where x which 
refers to the right-most site for the y-th "row" of sites 
from the bottom of the system block. By comparison, 
the index x represents the x-th left-most site for the y- 
th "row" of sites from the top of the environment block 
in this case. Thus, system and environment blocks with 
the same indices refer to block shapes which are related 
by a transformation of a rotation of 180° (about some 
arbitrary central point). 

The density matrix for the free spin SI and the sys- 
tem block is obtained by integrating out the degrees of 
freedom of the other spin S2 and the environment block. 
The new Hamiltonian and spin operators are then deter- 
mined in terms of this "augmented block" of the system 
block and its neighbouring spin 51. All of the informa- 
tion regarding these system blocks of shape (x,y) (for 



FIG. 4: The first step in the infinite-lattice approach. Note 
that the positions of the sites in the bottom-most row of the 
environment block are reversed from left to right. "Free" spins 
SI and S2 are again indicated by the unfilled circles. 



example, the number of states, the Hamiltonian, and the 
spin operators) is saved to disk for both of bottom and 
top blocks. This information is used later in order to 
start the full 2D lattice infinite-lattice algorithm. Hence, 
note that one firstly grows the width of the lattice until it 
reaches a preset size L and height 1, as shown in Fig. la. 
One then 'sweeps' from the left to the right such that the 
free spins S2 and SI run along the top of the ladder- like 
structure shown in Fig. lb. The quasi-2D blocks are now 
of height 2 and the process is stopped once the spin S2 
lies on the left-most site. 

In the first truly two-dimensional "step," the system 
block of index (x,y) = is copied into the environ- 

ment block, although the ordering of the sites on the bot- 
tom row of the environment block are reversed such that 
sites in the environment block run from L, L — l, 1, 
shown in Fig. ^for I2D = 2. The reason why one per- 
haps ought to do this is because it is hoped that it is 
beneficial to use a block of width L in order to retain 
the maximum number of bonds between the system and 
environment block (namely, L of them). Note that the 
system block is denoted by indices (3, 2) (again see Fig. 
E}. 

Presumably, the symmetries of this "reversed" block 
are slightly different from that of the unreversed case, 
and so this choice of environment block it is not ideal. 
However, the goal of this approach is to use previously 
defined blocks in order to quickly build up the final lat- 
tice in order for the finite-lattice algorithm to find the 
fixed-point of the RG process. However, note that that 
there may exist many alternate ways of performing this 
crucial "first step" and more is said of this below. Now 
that the first two-dimensional step has been described in 
detail it is possible to present the full algorithm. Note 
that we start from step 2 after defining the initial step 
mentioned above for the augmented block (2, 2) and that 
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FIG. 5: The second such "step" in the infinite-lattice ap- 
proach. Note new system block determined in the previous 
step, namely (2,?2d), is used in both the system and environ- 
ment blocks. 



the density matrix for the system block and spin SI and 
the relevant Hamiltonian and spin operators obtained for 
this augmented block are stored to disk at all points. The 
'steps' in the 2D algorithm is now given by: 

1. The system block of index (x,y) = (I^d — 

is now copied into the environment block and the 
ordering of the sites on the bottom row of the en- 
vironment block are reversed such that sites in the 
environment block run from L, L—l, 1, shown 
in Fig. 0] The two "free" spins are put at po- 
sitions (1,1 2 d) and (2, I^d) and the system block 
is defined by the left-most site in the second row, 
namely, block (x,y) = (3, 1 2 d)- The density matrix 
for the system block and the spin 5*1 is obtained 
and the effective Hamiltonian and spin operators 
determined. The results are then saved for both 
the top and bottom blocks of this shape of the aug- 
mented block, namely, (x,y) = (2,1 2 d). 

2. The next 2D "step" shown in Fig. 0uses system 
and environment blocks both indexed by (2, 1 2 d), 
although the sites in the bottom row of the envi- 
ronment block are not reversed in this case. Again, 
the density matrix for the system block and the 
spin SI is obtained and the effective Hamiltonian 
and spin operators determined. This information 
is saved to disk for the augmented block(s) given 
by 

3. In the next step the system block is referred to by 
the index (1,1 2 d) and the environment block by 
(3, ho), as is shown in Fig. The value of 1 2 d is 
now incremented. For example, after the routine 
which determines the initial ID blocks (see Fig. |3J) 
has been completed then 1 2 d is increment from a 
value of 2 to 3. 



FIG. 6: The third "step" in the infinite-lattice approach. 



4. The fourth type of "step" (or, more accurately 
stated, a set of steps) is in one in which one 
"sweeps" through values of lm = {1, 2, • • ■ , L— 3}, 
shown in Fig. . The density matrices are de- 
termined at each step and the information for the 
system block (x,y) = (1 2 d,Iid) for both the top 
bottom 2D blocks is saved. 

Clearly, the four types of different step described above 
may be used in an analogous way in order to sweep left- 
wards. (Note however that one starts from step 1 from 
now on). Furthermore, the final N = L x L lattice is 
formed by repeatedly sweeping right and then left, and 
so on. Thereafter the information stored to disk may be 
used in order to implement a finite-lattice algorithm in 
which one sweeps through all lattice sites repeatedly until 
a fixed point of the RG process is reached. It is further- 
more noted that prediction of the wave function is possi- 
ble for both the infinite-lattice and the finite-lattice algo- 
rithms, and that this can reduce the number of Lanczos 
iterations by up to an order of magnitude. 

It is hoped that this method might, in principle, al- 
low one to add in the degrees of freedom associated with 
the spin SI into new 2D blocks in a very controlled and 
efficient manner. Furthermore, it is hoped that the use 
of effective system and environment blocks at all points 
might to help in simulating the effects of the surround- 
ing lattice spins on these free spins - which might be 
especially pronounced for larger lattices. 

The infinite-algorithm presented here is, perhaps, one 
of the simplest 2D algorithms which utilises previously 
defined fully 2D blocks in order to build up the 2D lattice. 
The research presented in this article (and that also pre- 
sented in Ref. 0]) however conclusively proves that one 
may construct DMRG infinite-lattice algorithms which 
utilise previously defined 2D blocks at all points. Note 
however that more elegant solutions which also contain 
this property might be constructed using similar ideas, 
as mentioned above. For example, an environment block 
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Environment Block (1 1D +3,1 2D -1) 

1 2 .. \ D k+1 L+2 ► 

• L-l L 



1 2 •• \ D 



SI S2 



L+3 



environment, S2 



environment, S2 



(Al) 



\o+l L+2 \ D +3 



L-l L 



System Block (1 ID ,1 2D ) 



FIG. 7: The fourth and final set of "steps" infinite-lattice 
approach. The arrow shows the direction of the sweeping and 
note that 1\d runs over all indices {1, 2, • ■ ■ , L — 3}. 



of shape (1, lm — 2) might be used in step 1 of the 2D 
algorithm presented here. Another alternative is to use 
a line of spins (or partial line of spins), which would be 
necessary for step 1 only. Both of these alternatives to 
the 'step 1' used above would obviate the need to reverse 
the ordering of the lattice sites on the bottom on the 
environment block in this step. 

Note again that the Heisenberg Hamiltonian is given 
by Eq. and that we wish to solve the Schrodinger 
equation, H\^) = Egl 1 ^), where 1^} is the superblock 
ground-state wave function. The Hamiltonian may also 
be broken into distinct parts, given by 



H = H ; 
+ H, 



system ~f~ -^system, SI H - -^system, S2 



system , environment 



'environment 



It is noted that (by far) the largest amount of compu- 
tational time is spent evaluating iTsystem.enviromnentl*) 
during each Lanczos iteration. Indeed, the number of 
operations involved in determining this scales as LM 4 , 
where M is the number of states retained at each DMRG 
iteration. 

All system/environment interaction terms for the L 
distinct bonds between the system and environment 
blocks are determined and the gathered together before 
the Lanczos algorithm is invoked, and this information 
is then saved to local (RAM) memory. The number of 
operations involved in determining system/environment 
interactions is thus reduced by a factor of L, although 
the memory usage is commensurately increased by a fac- 
tor of M which puts a limit on the maximum number 
of the states M. Note that future calculations will dis- 
tribute this memory usage up by the use of parallel pro- 
cessing, such that each processor which on a given part of 
the whole system block/environment block interactions. 
This not only reduces the memory usage of each individ- 
ual machine but also considerably speeds up the DMRG 
algorithm. 

Furthermore, the data for the lower and upper two- 
dimensional blocks for indices (x, y) is also saved to disk. 
This is because the memory required to store all of the 
information for the spin operators scales with L 3 M 2 
which again very quickly becomes prohibitive. However, 
as this data is read from disk once per DMRG itera- 
tion, this approach is not too inefficient. The problem 
of system/environment block interactions appears to be 
present in all 2D implementations of the DMRG method 
in 2D and so an efficient implementation of these terms 
is crucial to a good DMRG code. 
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